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ABSTRACT 



The quasar PKS 1510-089 (z = 0.361) was observed with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes during high states 
. . 'in the optical and GeV bands, to search for very high energy (VHE, defined as £ > 0. 1 TeV) emission. VHE y-rays were detected with a 
^ ' statistical significance of 9.2 standard deviations in 15.8 hours of H.E.S.S. data taken during March and April 2009. A VHE integral flux of 
k> ' ^(0-15 TeV< E < l.OTeV) = (1.0 ± 0.2s,„ ± 0.2^,^) x 10"" cm^^s"' is measured. The best-fit power law to the VHE data has a photon index of 
. r = 5.4 ± 0.7s,a, ± 0.3sys. The GeV and optical light curves show pronounced variability during the period of H.E.S.S. observations. However, 
^ ■ there is insufficient evidence to claim statistically significant variability in the VHE data. Because of its relatively high redshift, the VHE flux 
^ " from PKS 1510-089 should suffer considerable attenuation in the intergalactic space due to the extragalactic background light (EBL). Hence, the 
measured y-ray spectrum is used to derive upper limits on the opacity due to EBL, which are found to be comparable with the previously derived 
limits from relatively-nearby BL Lac objects. Unlike typical VHE-detected blazars where the broadband spectrum is dominated by non-thermal 
radiation at all wavelengths, the quasar PKS 1510-089 has a bright thermal component in the optical to UV frequency band. Among all VHE 
detected blazars, PKS 1510-089 has the most luminous broad line region (BLR). The detection of VHE emission from this quasar indicates a low 
level of 7 - 7 absorption on the internal optical to UV photon field. 
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1. Introduction 

Blazars are a composite class of active galactic nuclei (AGN), 
consisting of BL Lacertae-type objects (BL Lacs) and flat- 
spectrum radio quasars (FSRQs). They are differentiated by the 
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presence (FSRQs) or the absence (BL Lacs) of strong emission 
lines in their spectra. The broadband spectra of blazars are dom- 
inated by no n-thermal emission, char acterized by rapid variabil- 
ity (see e.g.. IWagner & Witzellll995h in all frequency regimes, 
with high and v ariable polarization in the radio and optical fre- 
quency regimes jAUer et al.ll2003l: iMead et al.lll990h . More than 
three dozen blazars have been detected in VHE y-rays, the over- 
whelming majority of which belong to the BL Lac class. 

PKS 1510-089 is a F SRQ at a redshift of z = 0.361 
dBurbidee & Kinrnar] 19661). with highly polarized radio and 
optical emission dStockman et al.l Il984l) . At the milliarcsec- 
ond scale individual com ponents in the radio jet show appar- 
ent superluminal m otion (iHoman et al ] l2001h as high as 46c 
(iJorstad et al.ll2005h indicating a small inclination angle to the 
line of sight and high bulk Lorentz factors. Very long baseline in- 
terferometry (VLBI) observations of this highly polarized radio 
jet shows large misalignment be tween the milliarcs econd and 
the arcsecond-scale components jHoman et al. I I2OOI . This can 
be explained by the high bulk Lorentz factor in the jet, which 
makes a small jet bending appear much larger in the projected 
orientation seen by an observer 

The broadband spectrum of this source has a synchrotron 
component that peaks b etween millimeter and IR wavelengths. 
iMalkan & Moord (Il986h report broad e mission lines in the spec - 
trum of PKS 1510-089 (confirmed bv iTadhunter et al.l Il993l) . 
as well as a clear UV excess (the "blue bump") on top of the 
non-thermal continuum. The blue bump is attributed to ther- 
mal emission from the accretion disk. The high energy com- 
ponent in the spectrum extends from soft X-rays to GeV y- 
rays. PKS 1510-0 89 has been extensively monitored in X-rays 
(e.g.. lJorstad et al.ll2006.) a nd is known to be a br i ght y-ray emit- 
ter fro m the EGRET era (iHartman et"aD Il999l) . iKataoka et aP 
(l2008l) have shown that the quasi-simultaneous broadband spec- 
tral energy distribution of PKS 1510-089 can be well described 
by an external Compton model with seed photons from a dusty 
torus. In the soft X-ray band, Suzaku data su ggest a hardening 
in the spectrum, which iKataoka et al.l (l2008h propose could be 
due to either a small contribution from the synchrotron self- 
Compton component, or from bulk-Compton scattered radiation. 
Fermi-LAI measurements of the high energy (HE, defined as 
100MeV< E < 100 GeV) spectrum of PKS 1510-0 89 in diff^er- 
ent flux states during 2008-2009 were presented in lAbdo et all 
(I2010ah . The average HE spectrum (derived from the entire data 
set presented therein) is well described by a log-parabola model, 
dN/dE oc (E/Eq)-''-'^^''^'^"^«\ with the following best-fit values 
for the three free parameters - the spectral-shape parameters 
a = 2.23 + 0.02 andy6 = 0.09 + 0.01, and an integral photon flux 
above 100 MeV of (1.12 ± 0.03) x lO^^'cm-'s"', with the ref- 
erence energy, £0, fixed at 300 MeV. In the multiwavelength 
data presented by lAbdo et al ] (I2010ah no correlation between 
the flux variations in the HE and X-ray bands is found. The 
authors report a positive correlation between the HE and the 
optical band. They find evidence for a 13 day lag between the 
HE and optical 7?-band light curves (with the HE light curve 
leading). They argue that this behavior can be used to rule out 



Table 1. H.E.S.S. data and analysis results 



Data Set 


MJD 




Live time 


On'' 


Off 


y 


S'' 






n 


(hours) 


(counts) 


(counts) 


(counts) 




All data 


54910-54950 


22 


15.8 


823 


5750 


248 


9.2 


4 Tel. 


54910-54918 


23 


6.7 


447 


3172 


159 


8.2 



a change in beaming as the main driver for variability in the 
source. |Abdo et al. (2010a) make an estimate for the mass of 
the black hole in this source of 5.4 x 10^ Mq, using a model 
for the accretion disk temperature profile and the measured UV 
flux. They deduce a maximum isotropic y-ray luminosity of 
^ 2 X 10"*^ ergs"'. Because both these estimates are nearly an 
order of magnitude smaller than the values they obtain for more 
distant FSRQs, such as 3C 454.3 or PKS 1502+106, they argue 
that PKS 1510-089 could be an atypical FSRQ. At the same 



Notes. The first row gives all data taken in stereoscopic mode (i.e., > 2 
telescopes taking data). The second row gives the subset where all 4 
telescopes were operational. 

The time-weighted average zenith angle. 

Signal + background events around the source position. 

Background events from ofF-source region. 

The statistical significance in standard deviations. 



time, the high ratio between y-ray and synchrotron luminosities 
is a typical feature of an FSRQ. 

The luminous optical-UV photon fields (broad line emission 
and the blue bump) in FSRQs can cause substantial absorption 
of VHE photons by electron-positron pair p r oduction (see, e.g., 
'Do nea & Proflieroel l2003t iLiu & Bail 12005 iPoutanen & Stern I 
i2oioir If the VHE emitting region were to be embedded within 
the broad-line region (BLR), immersed in the reprocessed 
accretion-disk emission, the VHE y-rays may not escape the sys- 
tem. On the other hand, VHE emission has been reported from 
two FSRQs, viz. 3C 279 dAlbert et al.ll2008h and PKS 1222+216 
dAleksic et al.ll201 lb . 

PKS 1510-089 is a good candidate for a VHE emitting 
FSRQ because it is a bright GeV blazar with a jet that shows 
highly relativistic bulk motion (hence beaming effects should 
be strong). PKS 1510-089 was observed with H.E.S.S. (High 
Energy Stereoscopic System) in March-April 2009 to search 
for VHE emission, when it was reported to be flaring in the 
HE and optical domains. MWL data from this p erio d have been 
published elsewhere, e.g. jMarscher et aD(l2010l) and lAbdo et al.l 
(I2010ah . The H.E.S.S. data are presented followed by the rele- 
vant optical and Fermi-LAiT data, used for triggering H.E.S.S. 
observations and in the discussion section. 



2. Observations and analysis results 

2.1. H.E.S.S. 

Following reports in March 200 9 of flaring activity in 
PKS 1510-089 in the HE domain dp 'Am mandoeTaD 120091: 
iPuceUa et ai]|2009l: IVercellone et al.] |2009) as well as in the op- 
tical frequencies recorded with the ATOM telesc ope (an opti- 
cal flar ing st ate was also independe ntly reported bv lVillata et al.l 
d2009l) and |Larionov et al. (2009|)), it was observed with the 
H.E.S.S. array (.Hinton 2004 : Aharonian et al.|[2006b.) between 
MJD 54910 and MJD 54923. Subsequent observations were per- 
formed between MJD 54948 and MJD 54950 following more 
optical flaring and another HE flare that was reported by 
ICutini et aP (|200p). A total of 1 5.8 hours (corrected for dead 
time) of data passing quality cuts dAharonian et al ll2006bl) were 
obtained, with zenith angles between 14° and 41°, from all ob- 
servations. 

The data were analyz ed using the Model Analysis 
dde Naurois & RoUandl l2009l) . The relatively high redshift of 
the source implies that the VHE flux should be strongly at- 
tenuated due to the extr agalactic backg round light (EBL) in 
the optical to IR regime dNikishovlll962h . The amount of EBL 
extinction increases with the energy of y-ray photons, result- 
ing in a steepe ning of the spectrum mea sured in the VHE 
band (see, e.g., ISalamon & Steckerl Il998l) . Therefore, loose 



2 



H.E.S.S. Collaboration: Discovery of VHE y-rays from the quasar PKS 1510-089 



w180: 

|l6o|r 

"^140. 



H.E.S.S. - PKS 1510-089 

Exposure: 15.8 live hours 

Excess: 248.0 y, Significance: 9.2 o 




20 IVIar 
2009 



30 IVIar 
2009 



9 Apr 
2009 



19 Apr 
2009 



29 Apr 
2009 



(degO 



Fig. 1. Distribution of the squared angular distance of y-ray candidate 
events around the position of PKS 1510-089. The angular distribution 
(in terms of the square of the angular distance between the source 
position and the reconstructed arrival direction) of events around the 
position of PKS 1510-089 is shown. On-source (signal+background) 
events are shown as hatched histogram, whereas the off-source (back- 
ground) events are shown as points with error bars. The on-source re- 
gion is defined by a 6^ < 0.0125 deg^, shown by the vertical dotted line. 
The dashed curve is the PSF modeled as a Gaussian, and a constant 
scaled to the average background level. 



cuts (Ide Naurois & RoUandl l2009l) . that result in a lower en- 
ergy threshold, are used in the analysis. A total of 248 y-ray 
candidates were recorded from the source direction (Tablell] 
first row), which corresponds to a firm detection wit h a statis- 
tical significance of 9.2 cr (following the method of iLi & Mai 
Il983h . The distribution of the squared angular distance of events 
around the position of PKS 15 10-089, determine d using the "re- 
flected background" method (iBerge et al.ll2007l) . shows a clear 
excess in the source region as compared to the background re- 
gion (Fig.[T]). The observed excess is consistent with point-like 
y-ray emission from PKS 1510-089. The best-fit position of 
the VHE y-ray excess is at RA = 1542"5Z2 + 1^8stat + l-3sys 
(J2000), Dec. = -9°6'21'.'6 ± 26'.'5,smt ± 2QI!0,,, (nQ00). Th is is 
compatible with the optical position dAndrei et al.l l2009l) of 
PKS 1510-089, at a separation of 33'.'3 ±22'.'3, within the statis- 
tical and pointing errors of H.E.S.S. 

The VHE light curve (in 1 day bins) is shown in the top 
panel of Fig.|2l The best-fit integral flux (> 0.15TeV), ob- 
tained by fitting a constant-flux model to the light curve, is 
(8.0+ 1.6) X lO^'^cm-^s-i (^2= 20.4 for 11 degrees of free- 
dom). The X' test for variability, i.e. for a null hypothesis that 
the flux is constant, yields a p— value ^0.11 for this value of 
the test statistic. This is therefore insufficient evidence to 
claim statistically significant variability (at a 99% confidence 
level). Because a few bins in the light curve have low statistics 
(fewer than 10 counts in the source region and/or background 
region), the appropriate p-value was derived from simulations 
of a large number (~ 10^) of light curves. An average flux of 
8.0 X 10"'^ cm"-s"' was assumed, and the actual exposure times 
and instrument effective area in each bin were taken into ac- 
count. All bins could be thus considered in the test for variability. 

To get the best spectral reconstruction, an additional 
quality criterion was applied, only those data that were 
taken with the full array of 4 telescopes were accepted (see 
Table[Tl second row). A total of 6.7 hours of good quality 
4-telescope data were taken during all observations, yielding 
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Fig. 2. Multiwavelength light curves of PKS 1510-089 for the period 
between MJD 54910 to MJD 54951 in terms of integral fluxes in the 
respective bands along with the Icr eiTor bars. Panel (a) shows the VHE 
light curve (one day bins) from H.E.S.S. The horizontal line is obtained 
by fitting a constant-source model to the data. Panel (b) displays the HE 
light curve derived from Fermi-hAT data. The black open circles show 
the integrated fluxes in one day bins, whereas a finer binning of 4 hours 
is shown in magenta. The solid line is the average level, and the dashed 
line is the threshold level used for deriving a flare spectrum during the 
high state around MJD 54917. In panel (c) the i?-band optical fluxes 
measured with ATOM are shown. 



159 y-ray candidates from the source direction and a statis- 
tical significance of 8.2 cr. The ene rgy spectrum is d erived 
using a forward-folding technique (iPiron et al.l 1200 lb . The 
analysis threshold, ^thi 0.15 TeV, is given by the energy 
at which the effective area falls to 10% of its maximum 
value. For these observations the maximum of the measured 
differential rate is also at this energy. The likelihood maxi- 
mization for a power-law hypothesis, dN/dE - NoiE/E^))'^, 
in the energy range 0.15 TeV-l.OTeV, yields a spectral 
index of F = 5.4 + 0.7stat ± 0.3sys and a normalization con- 
stant of A^o = (1-1 ± 0.2stat + 0.2sys) X 10"'° cm^^s^'TeV"' 
at the decorrelation energy, £0 = 0.18 TeV (equivalent to 
a X' of 10.3 with 7 degrees of freedom). The spectral 
slope is steep compared to other VHE detected blazars, for 
exam ple, cf. F = 4.1 1 + 0.68stat for 3C279 (lAlbert et al.l 
I2OO8I) . The reconstructed H.E.S.S. spectrum is shown 
in Fig.|3] The integral flux, 7(0.15 TeV <£< 1.0 TeV) 
= (1.0 + 0.2stat ± 0.2sys) X 10"" cm"^s"' corresponds to 
K 3% of the Crab Nebula flux in the same band. Fitting a 
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Energy (TeV) 



Fig.3. VHE spectrum of PKS 1510-089 measured with the H.E.S.S. 
instrument, during March 2009. The solid line is the best-fit power law 
obtained using the forward-folding method (see text for details). The 
butterfly is the 68% confidence band, and the points with errorbars (1 cr 
statistical errors) are the energy flux. Arrows denote the 99% C.L. upper 
limits. 

broken-power-law or a power law with an exponential cutoff 
does not give a better fit. A spectrum obtained from the entire 
H.E.S.S. data (first row in Table[T]i is compatible with the 
spectrum derived from the high-quality 4-telescope data, within 
the statistical errors. The results have been cross-checked using 
a different a nalysis method b ased on multivariate analysis 
technique (see lOhm et alj|2009l and t he references therein), an d 
an independent calibration procedure (lAharonian et al ] l2006bl) . 

2.2. Fermi-LAT 

It is desirable to compare the VHE spectrum to the Fermi- 
LAT (lAtwood et a l. 2009) spectrum during the contemporane- 
ous period. While the Fermi-hPC^ team has published a spec- 
trum for this source from this period, the integration period used 
(~ 1 month) is much larger than the shortest flaring timescale 
seen by Fermi dAbdo et al. l l20T0al) . The intention here is to ob- 
tain a spectru m for the highest sta te within the FermZ-LAT flares 
(as defined in lAbdo et al. (l2010ah ) that is contemporaneous with 
HESS observations. Thus the Fermi-LAT data were analyzed 
using the publicly available Fermi Science Tool^ (v9r23pl- 
fssc-20111006) and the P7SOURCE_V6 instrument response 
functions. The light curve over a contemporaneous period as 
the H.E.S.S. observations is produced by a binned likelihood 
analysis retaining photons (the source class events) with ener- 
gies between 200 Me V and lOOGeV from a region of inter- 
est of a radius of 10° around the position of PKS 1510-089. 
All sources from the Fermi-LA ^I two-years point source cata- 
log (2FGL) dNolan et all 1201 2|) within an angular distance of 
15° of PKS 1510-089 were modeled simultaneously. Pass 7 
model^ of the Galactic and extragalactic backgrounds were 
used. For these two diffuse backgrounds, the normalizations are 
treated as free parameters in the likelihood analysis. As men- 
tioned before the Fermi-LAT spectrum of this source is bes t 
described by a curved log parabola model (lAbdo et al ] l20T0al) . 
Therefore, a log parabola model is used for the source in the 
likelihood analysis, using the gtlike tool, with the normaliza- 

' |http://fenni. gsfc.nasa.gov/ssc/data/analysis/documentation/Cicer one7] 
^ [http://fenni.gsfc.nasa.gOv/ssc/data/access/l at/Back:gro undModels.htmll 



tion and spectral parameters left free. The light curve derived 
between MJD 54909.5 and MJD 54951.5 is shown in Fig.[2l 
middle panel. The average (daily binned) integral flux between 
200MeV-100GeV is (1.26 ± 0.03) x lO^^cm-^g-i. Two flares 
are evident, one centered around MJD 54916, and the second 
centered around MJD 54948. Whereas the H.E.S.S. observations 
during the first HE flare had good coverage and data quality, the 
H.E.S.S. observations of the second HE flare commenced only 
after the HE flare had peaked. A finer binning of 4 hours was 
used around MJD 549 16 to precisely trace the development of 
the flare. 

To characterize the range of variations in the HE band, 
spectra were extracted over two different epochs from Fermi- 
LAT data. A spectrum for the long-term average state was de- 
rived from the first two years of Ferm/-LAT data and a sec- 
ond spectrum was extracted for the high state centered around 
MJD 54916. The region of interest, the sources modeled together 
in the binned likelihood analysis, and the parametrization of the 
source and the background were identical to that used for the 
light curve generation. The long-term average spectrum was de- 
rived from the data taken between MJD 54682-MJD 55412 (i.e., 
between 04.08.2008-04.08.2010) for energies above 200 MeV 
using a binned likeli hood analysis. This results in a total Test 
Statistic (TS, e.g., see lMattox et"ai][T996l) of 29440.8. The Hke- 
lihood analysis was applied in an iterative way, such that the 
parameters of the 2FGL sources within 15° of PKS 1510-089, 
and the spectral-shape parameters of PKS 15 10-089 are fixed in 
successive steps, leaving only the normalizations of the source 
and the two diffuse background components free to vary in 
the last step. The best-fit values for the log-parabola model 
(dN/dE = «o [E/Eor""'^ in[£/£o]) are - a normalization of «o = 
(1.41 + 0.02) X 10"^ cm^^s^'MeV"', and slope parameters of 
a = 2.21 + 0.03 and /3 = 0.083 + 0.011. The choice of the ref- 
erence energy, Eq, does not affect the spectral shape and was 
thus fixed at 260 MeV. These pa rameters are consis tent with 
the average spectrum presented in lAbdo et ak l (l20T0al) . A spec- 
trum in logarithm bins in energy (using a log-parabola source 
model) is also derived for the same period. The average and the 
energy-binned spectrum are shown in Fig.|4] left panel. A Fermi- 
LAT spectrum that is strictly simultaneous to the H.E.S.S. spec- 
trum (that is derived from a total exposure of 6.7 hours, with 
gaps in the observations, as well as varying live times for in- 
dividual exposures) cannot be constructed, because an integra- 
tion period of at least a few days of Fermi observations is re- 
quired to derive a meaningful spectrum. Moreover, the H.E.S.S. 
data used for spectral analysis were spread over a period of 8 
days, MJD54910-MJD54918, during which large flux varia- 
tions were seen at the Fermi-LAT energies. Because nearly half 
of the H.E.S.S. exposure was duiing the brightest phase of the 
GeV-flare centered around MJD 54916, the Fermi-LPiT data dur- 
ing this period were used to derive a spectrum representing the 
brightest HE flux state, contemporaneous to the H.E.S.S. obser- 
vations. Data taken between MJD 54914.8 to MJD 54917.5 were 
used, where the flux points in the Fermi-LPiT light curve (in 4- 
hour bins) were above a threshold flux of 2.0 x lO^^cm^^s"' 
(see Fig.|2|i. This threshold flux was chosen to ensure adequate 
statistics to construct spectral bins reaching at least lOGeV. The 
analysis of this data-set yielded a total TS= 2354.0, with the 
best-fit values for a log-parabola model given by a normalization 
of no = (10.5 + 0.7) x 10"' cm^^s^'MeV"', and slope parame- 
ters of = 1.81 + 0.13 andyS = 0.161 + 0.048 (reference energy 
fixed at 260 MeV). The corresponding energy-binned spectrum 
is shown in Fig. |4| (left panel). This is roughly an order of mag- 
nitude brighter than the average flux, though it should be noted 
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that the spectrum is biased towards the low-energy bins because 
of the lack of statistics at higher energies. 

2.3. ATOM 

The ATOM telescope (Autom atic Telescope for Optical 
Monitoring), iHauser et al.l (l2004l) . is a 1 m class optical tele- 
scope, located at the H.E.S.S. site. PKS 1510-089 is regularly 
observed with ATOM. From the onset of the H.E.S.S. observa- 
tion campaigns, the frequency of optical observations was in- 
creased to four observations in the /?-band (~ 640 nm) per night. 
The /?-band light curve, with 500 s exposure per point and not 
corrected for galactic extinction, is shown in Fig.|2] The opti- 
cal light curve shows clear and pronounced variability (a re- 
duced = 55.3, with 117 degrees of freedom from fitting a 
constant-source model, and the highest flux deviating more than 
20 cr away from the period average). A prominent flare centered 
around MJD 54917 was observed, with another flare occurring 
after MJD 54925 that could not be followed up with H.E.S.S. ob- 
servations because it happened around a full-Moon phase when 

H. E.S.S. does not operate. A third brightening was seen around 
MJD 54946. 

These observations show PKS 1510-089 in a relatively high 
state compared to its average optical brightness, measured with 
ATOM, between May 2007 and August 2009. For example, 
in the 7?-band, compared to the long-term average flux of 

I. 862 ± 0.004 mJy, the average during the H.E.S.S. observa- 
tions was 2.597 + 0.009 mJy and a maximum measured flux of 
5.724 + 0.150 mJy. 



3. Discussion 

The VHE y-ray flux of distant (z > 0.1) blazars suffer sig- 
nificant extinction due to pair-production on the diffus e UV- 
IR photon field in the intergalactic medium (.Nikisho\^ 119621 : 
iGould & Schrederlll966h . This phenomenon can be used to de- 
rive upper limits on the photon densit y of the UV-I R part of the 
EBL (e.g., Aharonian et al. 2006a; M azin & Rauel l2007). 

Here upper limits are derived by comparing the VHE spec- 
tra measured with H.E.S.S. to an extrapolation of the con- 
temporaneous HE spectral measurement s obtained with Fermi- 
LAT , following the procedures in, e.g., Georgano poulos et al.l 
(I2OIO), Aleksic_e tal. (201 1) and.Meye r et al. (2012). Modifying 
Equation 2 from lAleksic et al.l ( 1201 lb . to include the effect of 
systematic errors, the 95% confidence upper limit on the EBL, 
expressed in terms of the optical depth, Tmax, is given by 



Tmax(£') = In 



(1 - e) X (Fobs(£) - 1.64 X AFobs(£)) 



(1) 



where Fo\,s{E) is the measured flux, AFobsiE) is the correspond- 
ing 1 cr statistical error, FimiE) is the assumed unattenuated in- 
trinsic flux and s is the systematic error expressed as a fraction 
of the measured flux. For a steep spectrum this systematic error 
(i.e. the uncertainty in the absolute normalization of the flux) is 
~ 25% (i.e. s - 0.25). The systematic error is conservatively 
considered as a net bias in the flux measurements. It is the factor 
by which FobsiE) could overestimate the true flux. The estima- 
tion of Tniax in Equation[T] therefore involves a correction factor 
of (1 - e). The extrapolation of the HE spectrum measured with 
Fenni-LAI into the VHE regime is considered as an estimate 
for the unattenuated spectrum. The large variations in the HE 
light curve during the H.E.S.S. observing period require the need 



of contemporaneous spectral measurements. As already men- 
tioned, a strictly simultaneous Fen«/-LAT spectrum correspond- 
ing only to the short exposure of H.E.S.S. (column 2 of Table[Tl) 
cannot be obtained because of low statistics in the Ferm;-LAT 
data. Even the spectrum obtained during the high flux state (inte- 
grated over 2.7d, centered around MJD 54916) does not provide 
sufficient statistics at the high-energy end. Therefore, a conser- 
vative estimate resulting in the highest possible upper limit on 
the opacity is derived by re-normalizing the long-term average 
Fenni-LAI spectrum to the level of the high-flux state (dashed 
curve in Fig|4] left panel). This is justified because the long- 
term spectrum gives the best available description of the HE 
spectral shape up to the highest GeV energies, and changes in 
spectral shape during flares are smalfl This scaled-up spectrum 
is extrapolated to the VHE regime to estimate the unattenuated 
flux, Fi„t{E), used in Equation[T] It is assumed that the spec- 
trum smoothly extends from the HE to the VHE band, with- 
out spectral breaks or cutoffs. Any intrinsic spectral breaks re- 
sulting in lower fluxes at higher energies would yield a lower 
opacity from EBL attenuation. These assumptions along with 
the above mentioned considerations (the inclusion of systematic 
errors on the H.E.S.S. flux and the selection of the intrinsic spec- 
tral shape consistent with the brightest HE flux state) give firm 
upper limits to the EBL opacity. Other sources of opacity, e.g., 
absorption due to internal photon fields from the BLR or the 
dusty torus (Wagner et al. 1995), would decrease the contribu- 
tion from EBL extinction to the overall opacity and thus the pre- 
vious statement still holds. The EBL upper limits derived using 
Equation[T]are shown in Fig. |4] rig ht panel. The EBL limits from 
this distant FSRQ and the limits in lAharonian et al.l(l2006al) . that 
are derived from BL Lac type objects that are relatively nearby, 
are at a comparable level. 

Owing to their higher bolometric luminosities FSRQs are 
detected up to larger distances than BL Lacs. Furthermore, 
due to the bright emission lines in the spectrum of these ob- 
jects, which the BL Lacs lack by definition, the redshifts of 
FSRQs can be accurately measured. Thus the potential of putting 
strong constraints on the EBL extinction is very promising with 
FSRQs. However, in the case of PKS 1510-089, other factors 
currently offset this potential. The H.E.S.S. spectrum turns out 
to be steep as expected for a high redshift source. This results 
in a slightly higher systematic error of 25% on the VHE flux 
measurem ents (cf. ~ 20% for a Cr ab-nebula-like spectrum, in- 
dex ~ 2.4. lAharonian et al.l (l2006bh ). Furthermore, this sharply 
falling spectrum makes it difficult to collect good statistics at 
energies > 0.5 TeV and hence the statistical errors in this en- 
ergy range are also high. The combined effect of the large sta- 
tistical and systematic errors in the VHE spectrum (denoted by 
AFobs and s respectively in Equation[TJ, because of the steep 
VHE spectral slope of this source, limits the effectiveness of the 
method used above. Thus only weak EBL limits could be de- 
rived. Another aspect that limits the ability to put strong EBL 
constraints is the uncertainty in the intrinsic spectral shape. For 
the observations presented here there is insufficient evidence for 
claiming variability in the VHE band. However, because the 
unabsorbed HE band is highly variable (time scale of hours) 
it is difficult to obtain an accurate description of the intrinsic 
spectral shape (Fi„t in EquationlTJ due to the typically multi- 
day integration time required, hence making it difficult to put 
stricter EBL constraints. For example, if the estimated Fint at the 



3 lAbdo et alJ l l2010ah find a harder when brighter trend in the Fermi- 
LAT spectrum when the flux above 0.2 GeV is > 2.4 x 10^'cm^^s"'. 
However, the change in the photon index is small, < 0.2. 
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Fig. 4. HE and VHE spectra of PKS 1510-089 measured with Fermi-LAT and H.E.S.S., and upper limits on the EBL induced opacity. Left: 
The filled butterfly shows the H.E.S.S. spectrum as a 1 tr confidence band; the blue and red butterflies s hows the spectrum corre cted for EBL 
absorption, using the low model for the EBL from lGilmore et afl (120091) and the upper-limit EBL model from lAharonian et al.l ( l2006ah respectively. 
The triangles show the long-term average Ferm/-LAT spectrum, with error bars denoting 1 cr statistical errors; the solid curve is the log-parabola 
model obtained from a binned likelihood analysis of these data (see text for details). The open squares show the Fermi-LfiS spectrum derived 
from the integration period between MJD 54914.8 and MJD 54917.5 (error bars illustrate 1 a statistical errors). Arrows denote upper limits on 
the flux for those bins where insufficient statistics prohibits a flux measurement. The dashed curve is obtained by re-normalizing the long-term 
average Fenwi-LAT spectrum to the level of the high-flux state. Right: Upper limits on the EBL opacity derived from spectral measurements 
of PKS 1510-089 compared to published models. The solid line with arrows at both ends shows the maximum optical depth, corresponding to 
upper l imits on the EBL. For compari son the opacity corresponding t o the previously publi shed H.E.S.S. upper-limit EBL model lAharonian et al.l 
( l2006ah . the lFranceschini et al.l ( 1200 8h model and the low model from lGilmore et afl (120091) are shown. 



lowest energies were estimated to be 20% less than the value 
used, which is perfectly plausible when comparing to the fac- 
tor of ~ 4 variations in the Fermi fluxes seen during this period, 
the EBL UL from this source calculated as above wou ld have 
been more constraining than the lAharonian et af ] (l2006ah limits. 
With more multi-wavelength monitoring of bright FSRQs, such 
as PKS 1510-089 and other FSRQs that have less steep VHE 
spectrum, it could be possible to obtain a HE and VHE spectrum 
when the fluxes in neither band varies. This of course is bet- 
ter done with more sensitive Cherenkov telescopes, such as the 
H.E.S.S. II telescope array, which should provide richer statistics 
due to their higher sensitivities and larger energy coverage. This 
should allow putting stronger constraints on the EBL extinction 
by comparing their quiescent-state HE and VHE spectrum. 



The luminous broad line emission in the optical-UV band 
as well as the thermal UV excess in PKS 1510-089 indicate an 
intense internal photon field. This can in principle cause sub- 
stantial absorption of VHE y-rays within the BLR radius, or 
due to the reprocessed disk emission from the dusty torus at 
larger distances from the accretion disk. The H.E.S.S. detection 
of this source implies a low optical depth in the VHE emit- 
ting region. This can be tentatively explained by hypothesiz- 
ing that the VHE emitting region is in an optical ly thin part of 
thejet, presumably far outside the BLR, see e.g.. lTanaka et al. I 
( |20I Ih and Tavecchio et al. (2011) for discussion on the FSRQ 
PKS I222-H216. It should be noted that in the study made in 
iPoutanen & Stemi (l20I0h . where a search is made for signa- 
ture features of BLR absorption in the HE spectrum of several 
of the brightest FSRQs detected by Fermi, the HE spectrum of 
PKS 1510-089 is consistent with a negligible amount of absorp- 
tion due to line emission in BLR clouds. This is consistent with 
the hypothesis that the y-ray emitting zone is in an optically thin 



region. A detailed discussion of the internal opacity is beyond 
the scope of this work. 

4. Summary 

VHE emission was detected from the flat-spectrum radio quasar 
PKS 1510-089 with a statistical significance of 9.2 standard de- 
viations in 15.8 hours of H.E.S.S. data taken during March 
and April 2009. An integral flux, in the energy regime be- 
tween 0.15-1.0 TeV, of (1.0 + 0.2„at ± 0.2,y,) x lO"" cm^^g-i 
was measured, which is ^ 3% of the Crab Nebula flux. 
The spectrum is extremely steep with a photon index of 
r — 5.4 + 0.7stat ± 0.3sys for a power-law hypothesis. 

There is insufficient evidence to claim significant variabil- 
ity in the H.E.S.S. data. However, the multi-frequency data on 
PKS 1510-089 during the H.E.S.S. observation period shows 
clear and pronounced variability in the HE and optical bands. 
Using both the fermZ-LAT and H.E.S.S. y-ray spectrum, from 
200 MeV to 1 TeV, upper limits on the optical depth due to the 
EBL are derived. The EBL-Umits in this w ork are comparable 
with the limits in lAharonian et al. I (l2006al) . that were derived 
from BL Lac objects that are relatively nearby. FSRQs, such 
as PKS 1510-089, due to their higher luminosities compared to 
BL Lacs, can in principle, allow us to probe the EBL density to 
relatively higher redshifts. With the upcoming H.E.S.S. II tele- 
scope array more precise measurements of the spectra of distant 
blazars are expected, which will allow us to put stronger con- 
straints on the EBL. 
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